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We perform SPH+N-body cosmological simulations of massive disk galaxies, including a formalism 
for black hole seed formation and growth, and find that satellite galaxies containing supermassive black 
hole seeds are often stripped as they merge with the primary galaxy. These events naturally create a 
population of "wandering" black holes that are the remnants of stripped satellite cores; galaxies like 
the Milky Way may host 5 - 15 of these objects within their halos. The satellites that harbor black hole 
seeds are comparable to Local Group dwarf galaxies such as the Small and Large Magellanic Clouds; 
these galaxies are promising candidates to host nearby intermediate mass black holes. Provided that 
these wandering black holes retain a gaseous accretion disk from their host dwarf galaxy, they give 
a physical explanation for the origin and observed properties of some recently discovered off-nuclear 
ultraluminous X-ray sources such as HLX-1. 

Subject headings: galaxies : formation — galaxies : evolution — galaxies : halos — black hole physics 



1. INTRODUCTION 

Recent evidence for the existence of intermediate mass 
black holes (IMBHs) raises questions about how such 
objects might form and evolve. IMBH c andidates exist 
in globular cl usters (lUlvestad~et al. 2007) , nearby bulge - 
less galaxies (|Filippenko fc Hoi 120031: iBarth et "all 120041 ). 
and active galactic nuclei {Greene fc Hoi 120041 ) . Addi- 
tionally, off-nuclear ultraluminous X-ray sources (ULXs) 
have become increasingly promis i ng IM BH candidates 
(|Farrell et all [20091: Uonker et all [20Toh - A source of 
IMBHs may be the seeds of supermassive black holes 
(SMBHs) formed at high redshift; any seed that does 
not grow into a SMBH would today be observed as 
an IMBH. While the precise mechanism for seed black 
hole (BH) formation is unknown, there are several pos- 
tulated theories. One possible mechanism for SMBH 
seed formation is the direct collapse of pristine, low an- 
gular momentum gas, which forms BHs with mass on 
the order of 10 4 - 10 6 M ? (iKoushiappas et~all 120041: 
lLodato fc NatarajariliZOOl iBegelman et al.1 120061 ). An- 
other possibility is that the seeds are the remnants 
of Po pulation III stars, with masses around 10 2 — 10 3 
M Q (|Madau fc Reesll200ll IVolonteri et afl 120031) . Al- 
ternatively, the first nuclear star cl usters may collapse to 
form BHs of mass 1000-2000 M^ (jDevecchi fc Volonteril 
12001) . 

We perform SPH+V-body simulations of massive disk 
galaxies to explore the evolution of seed BHs in a cosmo- 
logical context. While previous works have addressed 
the growth of centra l SMBHs dDi Matteo et al.l 120081: 
lOkamoto etafl 120081 : iBooth fc Schavd 120091 ) There we 
specifically focus on those which do not end up in galaxy 
centers. We include BH growth by gas accretion and 
merging, as well as radiative feedback in order to form 



a fully self-consistent picture of BH growth and evolu- 
tion. We find that the tidal stripping of galaxies con- 
taining SMBH seeds leads to a population of "wander- 
ing" BHs within the larger galaxy halo. Such a phe- 
nomenon has been predicted in the case of idealized 
simulations of merging, galaxies (jGovernato et al.l 119941 : 
iKazantzidis et al.ll2005l) . though fully cosmological sim- 
ulations in the current ACDM paradigm are needed to 
test a suite of non-idealized scenarios. Wandering BHs 
may also be created by a grav itational slingshot from 3 - 
body black hole interactions ([Volonteri fc Pernal 1*2005) . 
gravitational recoil, or a distribution of Population III 
star remnants which have not spiralled into the central 
BH (|Schneider et al.ll2002D . In each of these cases, the 
times cale for dynamical friction is longer than a Hubble 
time (jTaffoni et al.ll2003l ). leaving a potentially substan- 
tial population of "wandering" IMBHs in galaxy halos. 

The recent discov ery of an off-nuclear IMBH candidate 
(|Farrell et al.l 120091 ) leads us to explore whether such an 
object can be explained by a "wandering" BH. The ob- 
ject HLX-1 is offset from the nucleus of its host spiral 
galaxy, and exhibits an X-ray luminosity of 10 42 ergs s _1 . 
This luminosity implies a lower limit to a black hole (BH) 
mass of 500 M , but the true mass may be much larger. 
We discuss under what circumstances HLX-1 would be 
observable. Preliminary results suggest that HLX-1 can- 
not be described by an isolated wandering BH, but a BH 
traveling with the remnant core of its parent galaxy may 
explain its observed properties. 

2. SIMULATIONS 
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We use the N-b ody code GASOLINE (IWadslev et al.l 
12004 lStadelll200 ID . an SPH tree code which incorporates 
star formation, gas cooling and hydrodynamics, and su- 
pernova feedback, and successfully models realistic galax- 
ies. Studies of cosmological simulations with GASOLINE 
have pro duced galaxies that follow the mass-metallicity 
relatio n dBrooks et aT"ll2007l) a nd the HI Tully-Fisher re- 
lation dGovernato et al.ll2009D . exhibit cold flow gas ac- 
cretion (jBrooks et al.l l2009) , and reproduce the distribu- 
tion of Damped Lyman Alpha systems at high redshift 
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TABLE 1 
Simulation Properties 



Run 


# within H v i r 




V circ 










(M ) 


(km/s) 


(kpc) 


(kpc) 


h239 


7408639 


8.32 x 10 11 


247.7 


242.9 


2.82 


h258 


7347383 


7.90 x 10 11 


225.4 


238.8 


3.81 


h277 


6624914 


7.09 x 10 11 


262.7 


230.3 


2.74 


h285 


7373714 


8.18 x 10 11 


238.2 


241.5 


2.24 



( Pontzen et al.l |2008|) . Additionally, iGovernato et al.l 
(2010) have shown how high resolution and physically 
motivated supernova feedback allow for the formation of 
a bulgeless dwarf galaxy with a dark matter core. The 
ability to form realistic disk galaxies is critical to our 
analysis of BH physics, because of the need to accurately 
trace the angular momentum of inspiraling gas which 
comprises the disk and eventually fuels the BH, as well 
as the star formation histories responsible for dispensing 
metals into the ISM via supernova feedback. 

For this Letter, we have simulated four different Milky 
Way-mass halos, which were selected from a uniform, 50 
Mpc volume and resimulated at a high r esolution using 
the v olume renormalization technique (Ka tz fc White! 
1993). All simulations are run with gravitational soft- 
ening e = 0.3 kpc, gas particle masses of 2.28 x 10 5 M Q , 
dark matter particle masses of 1.26 x 10 5 M Q , a Chabrier 
initial mass f unction dChabri cr 2003), and a WMAP year 
3 cosmology (|Spergel et al.ll2007f ) . Star for mation and su- 
perno va re cipes are described in d etail in iStinson et al.l 
(2006) and IGovernato et al.1 (|2007f) ; we adopt parameter 
values c = 0.1 and eSN = 1.0. Stars are allowed to form 
when gas reaches a threshold density of 1.0 amu cm -3 , 
which allows us the most physical representation of star 
formation given our resolution. While our dwarf galaxy 
satelli tes are not resolved to the level of IGovernato et aU 
•20101. for our purposes here this resolution is sufficient. 
Simulation properties are described in Table 1. Circular 
velocity V C i rc is measured using the wi dth at 20% of the 
peak of the simulated HI line profile fsee lGovernato et al.l 
2009). The disk scale length r s is determined by simul- 
taneously fitting an exponential + Sersic profile to the 
projected edge-on stellar surface density of the galactic 
disk. 

There is much uncertainty regarding the formation 
of the initial "seed" black holes which grow to become 
SMBHs; however, they must form ea rly on and grow 
quickly in order to form z ~ 6 quasars (|Fan et al] 1200 If ). 
We maintain a physically motivated scenario by assum- 
ing seed black holes form via t he direct collapse of ex - 
tremely low metallicity gas as in iBeeelman etaLl ([2006). 
However, our choice of seed formation scenario is not 
crucial. We focus here on the dynamical evolution of 
seed BHs, which will be the same regardless of where 
we set the initial mass (i.e. 100 M Q for the first stars 
or 10 5 Mq for direct collapse). The dynamical friction 
timescale is longer than the Hubble time even for these 
relatively heavy seeds, given the local densities involved 
here. Seeds are allowed to form if the pare nt gas particle 
meet s the criteria for star formation (see IStinson et al] 
l2006f) and additionally has zero metallicity. We designate 
a probability that a newly formed star will instead be- 



come a seed black hole with mass Mbh = 2.28 x 1O 5 M . 
This probablility is set to 0.10 in order to reproduce the 
predic ted black hole seed ha lo occupation probability at 
z = 3 (jVolonteri et al.ll2008f ) . Allowing seed BHs to form 
out of only pristine gas results in the truncation of seed 
BH formation at a redshift of ~ 3.5 due to the efficient 
diffusion of metals produced by the first supernovae into 
the ISM. This criterion causes seed BHs to form in the 
regions of early bursts of star formation, which are pri- 
marily the centers of the first massive halos to form in the 
simulation. In the event that more than one BH forms 
at the same time within a softening length, the BHs are 
merged into a single BH particle, with Mbh as the total 
of the individual masses. However, we make no assump- 
tions regarding the large-scale properties of the host halo, 
allowing the BHs to form based only on the properties 
of the local environment. BHs are not fixed at their halo 
centers, but are allowed to move as the forces upon them 
dictate. We minimize the effects of two-body interactions 
by adopting dark matter particle masses that are simi- 
lar to gas particle masses in the high resolution region, 
which greatly helps the central BHs stay at their halo 
centers. 

Black holes are presumed to accrete gas isotropically, 
following the formula for Bondi-Hoyle accretion. The 
accretion rate is limited by the Eddington rate, assum- 
ing a 10% radiative efficiency. In addition, a fraction 
of the rest mass energy of the accreted gas is converted 
to radiative energy, which i s then isotropically imp arted 
onto the surrounding gas (|Di Matteo et al.l 12005). To 
dissipate the feedback energy in a realistic way, we use 
a blast-wave feedback appr oach similar to the su pernova 
feedback recipe described in lStinson et al] (|2006f ) . We set 
the fraction of radiated energy given to the surrounding 
gas to be 0.1%. 

Black holes are allowed to merge if they (a) are within 
one another's softening length and (b) fulfill the criterion 
\Av 2 < Aa- Ar, where Av and Aa are the differences in 
velocity and acceleration of the two BHs, and Ar is the 
distance between them. 

We identify galaxies and their halos with AHF 
(|Knebe et al.l [200lt IG111 et al]|200l . which identifies a 
virial rad ius based on the overdensity criterion for a flat 
universe (|Grossl[T997f ). For each output, we identify ev- 
ery halo and all of the particles it contains; thus, we can 
track a BH's halo residency at each timestep. We use this 
process to define whether a black hole is stripped from 
its parent halo to reside in another. A BH is denoted as 
"stripped" if its parent halo overdensity can no longer be 
detected as a separate halo. The BH is then identified as 
a resident of the larger galaxy, most likely in the outer 
regions. 

3. RESULTS 

We provide a dynamical mechanism to place massive 
BHs in galaxy halos, which is simply a consequence of 
hierarchical merging in the ACDM paradigm. Detect- 
ing these objects may help constrain the mechanism of 
seed BH formation, and we investigate a few scenarios in 
which a wandering BH may be observed as a ULX. 

3.1. Wandering Black Holes 

For our four simulated halos, we show the radial distri- 
bution of BHs at z = 0, defined as the distance between 



Wandering BHs in Disk Galaxy Halos 



3 



the halo center and the BH, in Figured] Each of the four 
galaxies has a central BH and between 5 and 15 "wan- 
dering" BHs whose distances range from 10 - 100 kpc. 
The vast majority of wandering BHs have grown in mass 
by less than 2% since their formation. A small num- 
ber (5 out of 36 total) have undergone BH — BH mergers 
early in their lifetimes and thus have larger masses, but 
only one of these has experienced any substantial accre- 
tion events. In all cases accretion is effectively quenched 
when the host galaxies are torn apart and the BHs are 
left to wander throughout the halo, and thus any BHs 
present here are at or within a few factors of their origi- 
nal seed mass. 
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Fig. 1. — Top Panel: Distribution of black hole radial distances 
to their halo centers for 40 black holes in four simulated galaxies. 



Local Group d warf galaxies are promising t argets for 
IMBH searches (|Van Wassenhove et aLl l2010f l. In Fig- 
ure [H we show the masses of BH-hosting satellite galax- 
ies before the stripping process has begun. This distri- 
bution peaks between 10 9 — 10 10 M Q ; at halo masses 
of 3 x 10 8 M & the satellites become increasingly dark. 
Since very few stars form in these halos, it is unlikely 
that we are underestimating the number hosting mas- 
sive BH seeds, given our scenario. We estimate the ab- 
solute magnitudes of the sa tellites in these mass ranges 
using the STARBURST99 (jLeitherer et al.lll999D code, 
using stellar ages and metallicities from our simulations. 
This mass range corresponds to galaxies with magnitudes 
fainter than -15 in the V band, which includes several Lo- 
cal Group dwarfs, including the Large and Small Magel- 
lanic Clouds. Examining these local dwarfs for IMBHs 
may help constrain the locations and masses of SMBH 
seeds. A confirmed detection of an IMBH would pro- 
vide an upper limit to the initial mass of SMBH seeds 
and possibly allow us to differentiate between the various 
proposed formation mechanisms of such seeds. 

3.2. Connection with Off-Nuclear ULXs? 

While wandering BHs may be present in the Milky 
Way halo, observing such objects would require either a 
fortuitous gravitational lensing detection, or a triggered 
accretion event, resulting in X-ray emission. We inves- 
tigate situations where a wandering BH could undergo 
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Fig. 2. — The mass distribution of BH-hosting satellites before 
they are stripped by the primary. The solid line indicates the 
total mass of the satellites, while the dashed line is the stellar 
mass only. Our simulations are incomplete for halos below a total 
mass of 3 X 10 8 Mq, however due to the UV background field 
IjHaard t & Madau 1996:) most halos with mass below 10 s Mq have 
very few stars. 



substantial accretion in a galactic halo, which could ex- 
plain the existence of objects such as HLX-1. 
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Fig. 3. — Each curve represents the radial distance vs. time of 
a black hole to the halo center for the simulated galaxy h239. The 
red dashed horizontal line denotes a distance of three scale lengths 
(r s = 2.8 kpc) from the center; each instance of an orbit crossing 
this line is defined as a disk passage. The black dashed horizontal 
line shows the gravitational softening, while the dotted line repre- 
sents the virial radius. The yellow line indicates the central BH, 
while the green line denotes an inspiraling BH which has not yet 
reached the halo center. The other curves represent BHs which 
reside predominantly in the halo. 



While the maj ority of ULXs are a ttributed to star for- 
mation regions (jSwartz et al.l 120041) , a smaller number 
of off-planar ULXs are not cospatial with star forma- 
tion, and we focus on the possible origins of these objects 
here. We first examine whether a wandering BH passing 
through a dense portion of the host galaxy disk could re- 
produce the properties of off-nuclear ULXs. We trace the 
orbits of each BH and examine the frequency with which 
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they cross the disk region, which we define as three disk 
scale lengths (r s ). The last few Gyr of orbital history 
for the BHs in our simulated galaxy h239 are shown in 
Figure El Several of the BHs pass through the 3r s limit, 
which we denote as a disk passage. The average BH disk 
passage rate for all of the simulations combined is 10.6 
per Gyr; however, this rate increases with time (13.3 per 
Gyr for the last few billion years of galaxy evolution) due 
to the increased population of wandering BHs in the halo 
(since more mergers have occurred). 

The likelihood of observing such a disk passage de- 
pends on the column density of gas in the disk, the num- 
ber of passages made, and the velocity of the BHs. We 
use results from HI me asurements of the surface densities 
of local disk galaxies (|Lerov et al.l 120081 ) to estimate an 
average disk column density of 10 M Q pc -2 . The ability 
of the passing BH to attract a sufficient amount of gas 
to be observed as a ULX depends most strongly on its 
velocity, which ranges between 154 - 662 km s _1 , and on 
average is ~ 470 km s _1 for BHs at pericenter. To esti- 
mate a BH's luminosity, we developed a simple scenario 
where the passing BH is able to collect the gas within a 
radius of influence determined by its mass and pericenter 
velocity: M coU = Z disk Tr(GM B H /v 2 BH ) 2 . We perform a 
Monte Carlo simulation taking into account the observed 
ranges of BH velocities, and estimate that the amount of 
mass collected by a BH during its disk passage ranges 
from 10 -3 — 10~ 4 Mq. Assuming an accretion luminos- 
ity L ~ rjMc 2 with 77 = 0.1 and the canonical luminosity 
of Ljjlx — 10 39 ergs s _1 , the BH will accrete and ra- 
diate in a high state for a mean duration of 2000 years, 
though a slower BH may radiate for up to 10 5 years. 
However, if the ULX transitions from a high state to a 
low state at any point, the accretion rate would drasti- 
cally decrease and the observed emission could continue 
for much longer. 

This simple scenario, however, is not sufficient to re- 
produce the properties of HLX-1, which is at least 1 kpc 
away from the plane of its host galaxy. A BH passing 
through the disk moving at a few hundred km s _1 will 
only travel a distance of ~ tens of parsecs during the 
duration of the predicted accretion event. Such an event 
would appear to be a disk ULX source, possibly indis- 
tinguishable from those cospatial with star formation re- 
gions. HLX-1 is also 1000 times more luminous and 100 
times less massive (Wiersema et al. 2010., in prep.) than 
the assumed values in our calculation above. Taking all 
of these factors into account, it is unlikely that an iso- 
lated wandering BH passing through a galaxy disk can 
reproduce the properties of off-nuclear ULXs. 

However, our wandering BH scenario would be feasible 
if HLX-1 is not an isolated BH, but is traveling with a 
bound clump of gas and stars. The magnitude and colors 
of the detected optical cou nterpart of HLX-1 are consis- 
tent with a globular cluster ()Soria et a l. 2009). This clus- 
ter may in fact be the stellar remnants of a stripped dwarf 
galaxy core that is still bound to the BH. HLX-1 could be 
similar to the object Gl in M31, a globular cluster that 
may actually be the nucleus of a stripped dwarf galaxy 
(jMevlan et al.|[200ll ) and harbors t he most promising Lo - 
cal Group candidate for an IMBH (jUlvestad et al.ll2007t ). 
If we presume that the wandering BHs in our simulations 
retain a nuclear star cluster and gas reservoir from their 



parent halos, then the instance of a wandering BH pass- 
ing near the center of the primary could cause instabili- 
ties in its accretion disk, triggering an accretion event of 
sufficient magnitude to power a ULX with a high lumi- 
nosity as seen in HLX-1. Our simulations do not have 
sufficient resolution to follow the tidally stripped cores of 
galaxies in detail, though other simulations have shown 
that a tidally stripped dwarf galaxy ca n retain its core af- 
ter a close passage with the primary (jMaver k, Wadslevl 
2004). Thus, in the likely instance that a wandering 
BH retains the core of its host galaxy, its passage near 
the galaxy disk can explain the origin and properties of 
HLX-1. 

Previous studies have estimated luminosities for mas- 
sive BHs wandering through the ISM, but prior to this 
Letter n one have exp lored the issue in a cosmological 
context. iKrolild (|2004t ) showed that IMBHs with masses 
ranging from 10 2 — 10 4 M can produce luminosities of 
ULXs if they pass through or near molecular clouds. 
SMBH gravitational recoil events passing through the 
disk may exh ib it X-ray emiss i on of L > 10 39 ergs s 
(jFuiital l2009h . iMapelli et all pOOl performed an N- 
body+SPH simulation of IMBHs in a galaxy merger, and 
found that a few halo IMBHs reside in orbits that pass 
through the disk, which in our scenario may be observ- 
able as ULXs. 

4. SUMMARY 

We provide compelling dynamical scenario for the pres- 
ence of "wandering" massive BHs in the halos of galax- 
ies. A natural consequence of the hierarchical build up of 
galaxies in a ACDM scenario, the tidal stripping of galax- 
ies containing seed BHs can populate the halo of a mas- 
sive disk galaxy with wandering BHs. These objects of- 
ten retain their original seed mass, are found throughout 
the galaxy halo, and may pass through the galactic disk 
at an average rate of 10.6 Gyr" 1 . We predict that Lo- 
cal Group dwarf galaxies such as the Magellanic Clouds 
are likely to host IMBHs. Detections of these wander- 
ing BHs may give an upper limit to the initial mass of 
BH seeds, and may allow us to differentiate between the 
various proposed formation mechanisms of such seeds. 
Our scenario provides a physically motivated explana- 
tion for off-nuclcar ULXs as IMBHs which have been 
stripped from their host galaxies, if they retain a gas 
reservoir/accretion disk that, when dynamically destabi- 
lized, might be funneled toward the black hole and form 
an accretion disk. 
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